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In the framework of SU(2) singlet down type vectorlike quark model, we present a comprehensive 
analysis for decays B — > X a j, B — > X s l + l~ and B — > K-k. As for B — > X s j, we include the 
QCD running from the mass of the down-type vector quark D to weak scale in the scenario with 
the D quark much heavier than weak scale, and find that the running effect is small. Using the 
recent measurements of B — > X s l + l~, we extract rather stringent constraints on the size and CP 
violating phase of z s b, i.e., the tree level FCNC coupling for b — » sZ. Within the bounds, we 
investigate various observables such as forward-backward asymmetry of b — > sl + l~ , the decay rates 
of B — > X a "/, and B — » K-k. We find that (1) The forward-backward asymmetry may have large 
' derivation from that of the SM and is very sensitive to z s t,, and thus can be useful in probing the 

new physics. (2) By taking experimental errors at 2a level, both experimental measurements for 
B — * X s l + l" and B — > Kk decays can be explained in this model. 
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I. INTRODUCTION 

The study of flavor changing neutral currents (FCNC) in particle physics phenomenology has played a key role 



(N 
> 

in advance of high energy physics in the past decades. Due to the GIM mechanism, FCNC in the standard model 
(SM) arises only at higher loop level, and thus, it makes FCNC phenomena a privileged ground to search for signs 
of new physics beyond the SM. However, in extension of the SM such as vector quark model (VQM), the CKM 
• matrix is necessarily non-unitarity, leading to interaction Zsb at tree level, and hence potentially large new physics 
' contributions can be expected. 

The rare radiative decays B — > X s j and B — > X s l + l~ are sensitive probes of new physics [1]. The branching ratio 
O ' of the radiative decay B — > X s j has been measured by BaBar [2], CLEO [3], and ALEPH [4] and is in good agreement 
with the SM predictions [5,6]. Recently, the rare decays B — > X s l + l~ (I = e,fi) also have been measured by BaBar 
[7] and Belle [8] . The average value is [9] 



Br ex (B -> X s l + l~) = (6.2 ± l.l±i j) x 1(T 6 . (1) 



KJi, \ Although the process B — > X s ^ constrains the parameters of the VQM [10,11], since vectorlike down-type quark 
contributions to b — > 57 just occur at loop level as the case of the SM, the constraints on z s t,, the tree level FCNC 
coupling for b — > sZ, from B — > X s j are less restrictive compared to those from processes governed by b — > sl + l~ 
transition. 

There are some studies regarding the constraints on model with extra singlet quark [10-13]. In light of the improve- 
ment in the experimental data, it is necessary to present a comprehensive analysis in this model. Also, from the point 
of view of the model builder, it is important how the presence of the singlet quark may have impact on low energy 
phenomenology. In particular, mass of the singlet quark, the coupling to ordinary quarks and weak gauge bosons are 
very important issues. We extend the previous studies and take the following points of the VQM into account: 

(i) In the previous studies [10-14], the down-type vector quark D is integrated out with W, Z bosons and top quark 
together at mw scale, neglecting the QCD running from to weak scale. In this work, we also consider the scenario 
with the D quark much heavier than weak scale. Firstly, the down-type vector quark is integrated out, generating an 
effective six-quark theory at scale. By using the renormalization group equation (RGE), the effective field theory 
is run down to the weak scale, at which the W, Z bosons, Higgs and top quark are removed. Finally, the effective 
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field theory is running down to the b quark scale, as usual done in the SM. From viewpoint of the theory, if there are 
different scales in a model, including QCD running from heavier scale down to lighter one is important. 

(ii) For inclusive decay B — > X s l + l~ , we include the four-quark operator contributions to one-loop matrix elements 
of operator C 9 due to the tree level Zsb interaction. We also consider the long distance contribution from resonance 
4>. This is because the tree level FCNC b — ► sss generates the new decay chain B — > X s (f> — > X s l + l~ . We show that 
the dilepton mass distribution can be affected. 

(iii) With the various improved treatments in both short and long-distance contributions at hand, we obtain the 
rather stringent bound on the tree level Z FCNC coupling and its CP violating phase. Within this bound, we study 
various observables such as forward-backward asymmetry of b — > sl + l~ , the decay rate of B — > X s j. 

(iv) . The large electroweak penguin contribution to b — > s transition has been suggested in the present data of 
B — > i^7r [15-19]. In this work, we use factorization approach and study whether the large electroweak penguin 
contribution in B — > if-^ can be explained or not within FCNC constrained by rare decays b — > sl + l~. 

This paper is organized as follows: In Section II we give a brief description of the VQM. We present calculation of 
b — > S7, 6 — > sl + l~ transition in the VQM, including QCD running from down-type vector quark mass scale to weak 
scale in Section III. Some new operators are introduced and the contributions of electroweak penguin operators are 
taken into account. In Section IV , we extract constraints on size and phase of z s b from B — ► X s l + l~ experimental 
measurements. Using the bounds in Section IV, we evaluate the forward-backward (FB) asymmetry and its zero-point 
So for the rare B dileptonic decay, as well as the branching ratio of B — > X s -f, and show how they are affected by the 
new physics. Section VI contributes to the study of B — > Ktt decays. We found both experimental measurements for 
B — ► X s l + l~ and B — > Ktt decays can be explained within the framework of vector-like quark model. The anomalous 
dimension matrices needed in solving Wilson coefficients and all loop functions are collected in Appendix A and 
Appendix B, respectively. 



II. VECTOR QUARK MODEL 



In this section, we summarize the parameterization of quark mixings in vectorlikc quark model. We focus on the 
model including a singlet vectorlike down-type quark denoted by D, added to the standard model. 

The difference between the new quark and ordinary quarks of the three SM generations is that, unlike the latter 
ones, both left- and right-handed components of the former quark is SU(2) singlet. Then the CKM matrix Vckm is 
enlarged to 3 x 4 and can be expressed as [10] 

V CKM = ^2U ji V 0i *, (2) 

where U, V are 3x3 and 4x4 unitary matrices which relate the weak-eigenstates q^ to mass-eigenstates ql (q = u,d), 

4 = Vf> a dZ, ui = U^ul (3) 

The 4x4 matrix V covers ordinary (a — 1,2,3) and and vectorlike quark (a = 4). 

The fact that the vector like quark is isosinglet, leads to non-unitarity of the mixing matrix V as 

z af ^ ymySi " = E V ckmV^ km = S aP - V ai V fM *. (4) 

i=l i=l 

Geometrically, Eq.(4) for % = a ^ f3 = j shows the quadrangle in the complex plane [20]. The deviations from the 
standard unitary triangles (zij ^ 0) are going to vanish as the down type singlet mass (M) increases compared with 
electroweak breaking scale v. 

The interaction Lagrangian for the quarks with the W ± and Goldstone bosons \ reads 

LcC = V2W V ° KM "* {mu3L ~ md " R) dPx+ + ^ y cU/^^ + h.c, (5) 

while the interaction Lagrangian for the down- type quarks with the Z, Higgs H° and Goldstone bosons x° is given 
by [10] 
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C-NC = - 



cos 9 



w 



aJ3 



^z a/3 L + e d sin 2 9 w 8 a p 



d p Z u 



2m\y 



ct,0 



+i7T— Y\ d a z af3 [m da L ~ m d ?R} d x°, 



(6) 



where e d = — 1/3 is the electric charge of the down- type quark. 

So far, the results are general. To discuss the effects of singlet quark on unitarity of the CKM matrix, we must 
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keep the length of the sides of quadrangles to the order of j^. Then we can discuss the detailed structure of the 
quadrangles. For that purpose, we devise the parameterization of 3 x 4 non-unitarity matrix based upon the systematic 
expansion of jj. Using the parameterization, we show the quadrangle of bs sector, which will be used in later section. 
It is completely general to start with mass matrix of the down type quarks as follows [20] : 



M d 



with charged current 



mo J 
M 



u L -ff,K°d L , 



(7) 



where mo is a real diagonal 3x3 matrix with (m )ij = m iSij. J is a 3 x 1 matrix, J T — (Ji, J 2 , J3), and Ji is real 
while J 2 and J 3 are complex. M is a singlet quark mass parameter and can be taken as real. K° is a 3 x 3 unitary 
matrix and can be parameterized by A, A, p, r\ as in the CKM matrix of the SM. Note that due to non-vanishing 
Z FCNC couplings, the values of A, A, p, r/ can be different from those of the SM. The matrix A4 d Ai a can be 
diagonalized by a 4 x 4 unitary matrix W , 



W = 



n r 

S T 



WM d Mtw 1 = 



m 






M 2 



(8) 



Thus, the 3x4 non-unitary matrix Vckm is a submatrix of W, Vckm = (Q,R)- The matrices and S satisfy 
following equations: 



(mo + JJ r )n + MJS = Vim 2 , 
MJ^n + M 2 S = SM 2 , 
+ = 1. 



(9) 



By eliminating S and using R ~ -jj , we obtain £1 to the order of 



jj t 



M' 2 



± 2 
a 2 i 



as 

m 03 

^12 a "2" ^13 A 'i 



Am 021 



l 02 ^ 

A 



23 Am 2 



where 



J.J* 

' ^ m 0tj = m 0i 



in 



oj- 



(10) 



(11) 



Now we can contact the expression of Z coupling Zij with Ay in this approximation. From Eq.(4) and (10), we obtain 



Finally, considering m i ~ m d <C m 2 — m s <C m 3 ~ mj,, we have 



1 



Vckm — K° 



An 
2 



Z ds 







1 - 




(12) 



(13) 



where Si — arg(Ji) (i — 2, 3). We note that there are three independent CP violating phases. Thus, by parameterizing 
the CKM matrix, we link the CKM matrix Vckm to a unitarity matrix K° and Z FCNC couplings. 
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Using the experimental measurements for BdBd, KK mixings, CP asymmetry of B — > ^Ks and CKM matrix 
elements, we can constrain the Z FCNC in bd, sd sectors, and investigate the correlations among the Z FCNC in 
bd, sd, bs sectors. The detailed study of them will be presented elsewhere. In this work, we will focus attention on 
the Z FCNC in bs sector, assuming that the Z FCNC effects on mixings and decays of K, Bd are negligible. In this 
approximations, we obtained 

K ub K* us =A\\p-i V ), 
K cb K* s = AA 2 (l-y) + z sb , 

K tb K* ts = -A\ 2 {1 - ^) - AXHp - ill), (14) 
where 2£y = {VcK\i)ij f° r hj — lj 2, 3. The corresponding quadrangle in bs sector is shown in Fig. 1. 




FIG. 1. The quadrangle in bs sector. 



III. B — > 57 AND B -» St + r TRANSITIONS IN VQM 



In VQM, the down-type vector quark may be much heavier than weak scale. In a theory with different mass scales, 
the heavier scale should be integrated out firstly, then Wilson coefficients are run from heavier scale to low scale by 
using renormalization group equation. Only in case of mo is about the weak scale, can W, Z boson, Higgs boson and 
top quark be integrated out together. In this work, we consider two possibilities as follows: 



A. Scenario A: 8d = ml ^- H < 1 



In this scenario, we first integrate out the heavy D quark, introducing dimcnsion-5 and dimension-6 operators. By 
keeping only leading order terms of Sd, we obtain the effective Hamiltonian for b — > sj*(g*) as: 

W$?(6-7*Gf)) = ^-zlz^CiMOi. (15) 

* i 

A complete basis for the local operators is listed below: 

°LR = -TTT^mbSLTPbR, 
lD7T z 

0\ R = j^p e/2 g s m b s L a^T a b R G^, 
3 LR = Y^^^eedmhSLcr^bnF^, 
Qlr. = lv £ g 2 s m b H°H s L b R , 
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Qlr = ^ e 9 2 s m bX°X°SLbB,, 



16tt 2 
1 

2' 
.1 



R l H = W e g 2 a H H s L pb L , 



R 2 L H = i^g 2 s (d a H )H a s Lla b L , 



R 2 * = iH e gl(d°x )x SLi*b L , (16) 



where T a stand for SU(3) co i or generators, F^ v and G^ v are field strength of photon and gluon respectively. The 
covariant derivative is given by 

X> M = 0„ - i^ /2 g s G^T a - t^ /2 ee d A^ (17) 
with e = 4 - d. The tensor T£ vg (A = 1, 2, 3, 4) appearing in P L ' have the following Lorentz structures [21,22]: 

Tuva ~ dfivla-j 'Fui/a = 9p,cr^/vi F^ VIJ = Qva^lin = — * e /ifcrr7 75- (1$) 




7,9 




7,9 



FIG. 2. Matching conditions at scale mo in full theory (left) and in the intermediated effective field theory (right). <j>° can 
be H° and X °- 



To match the full theory into an effective theory, the diagrams to be matched at scale are displayed in Fig. 2. 
To determine the coefficients of the r bs<w \ T b ^^, r bs ^ (<£ = H°, X °) at scale m D , we just need to match the full 
theory with effective theory at tree level. We obtain 



C Q H R (m D ) = -C Q x R {m D ) = -jz 

y s 

i 

1 



C R i H {m D ) = 2C R 2H(m D ) = - — 



C R i, x (m D ) — 2C R 2, x (m D ) 



(19) 
(20) 
(21) 



Other coefficients can be determined by matching one loop diagrams shown in Fig. 2. After straightforward 
calculations, we have 
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C oiR (m D ) = (1 = 1,2,3), 
C p i,i(m D ) = C p i,3(m£)) = - 



11 

18' 



C p i,2(m D ) 
C p iA{m D ) 



9' 



C 2 PL (m D ) = 0. 



(22) 



The values are understood as sum of H°, x° contributions. Cancellation between Higgs and would-like Goldstone 
boson x° leads to Cq (trd) =0 (i = 1, 2, 3) in leading order. 

The running of Wilson coefficients Ci(fi) from mr> down to weak scale is governed by anomalous dimension 7 
through RGE 



(23) 



We calculate one-loop diagrams with operator insertions and present the anomalous dimensions 7 in Appendix A. 
Using the anomalous dimensions (Al) and (A2), we can solve the RGE (23) and have the coefficients of the operators 
at weak scale m^. They are given by, 



c o 2 LR ( m w) 
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247 



1096 
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1096 



336 . 113 

r 126 — r; 

8905 s 780 s 
168 ,113 

\ -Q 126 



8905 



168 ,113 

C 126 

8905 s 
247 



551 

223 ,8 1^2 

C 21 J ^3 

780 s 24 s ' 

223 5,2 3 ,16 
r 21 J m — —C21 

780 s 12 s 8 s ' 



C p i,i(m+ ) = C p i, 3 (m+ ) = -— C 21 - 



548 



791 
4932 







875 ,113 



2466 



C p i,4(m^) 
Cp 2 (m+ ) = 0, 



247, _jl 1 ,jl 14 ,113 
r 21 — — r 21 J r 126 

598 s 12 s 411 s : 



(24) 



where C = a s (m D ) /a s (m^). 

In order to continue running the basis operator coefficients from m\y scale down to b quark scale, we use the effective 
QCD-corrected Hamiltonian obtained by integrating out the W, Z bosons, would-be Goldstone boson, Higgs boson 
and top quark. The effective Hamiltonian describing b — > sl + l~ transition reads [5,23] 



H eff (b^si+n 

where 



4G F 
'^2 



KtbKL 
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Ci{n)Qi + C 77 ( A1 )Q 77 + CsgMQsg + C 9 (^)0 9 + C w (n)O w 

Qi = {siYLc^Cj^i^Lbi), 
Q 2 = (s^Lc)(c^Lb), 

Qs = (srLb)J2(Q^Lq), 
q 

Qi = {si^Lb^Y^iQjluLqi), 
<? 

Q 5 = (s^Lb)J2(q 7tl Rq), 

<i 

Q 6 = {sa^Lb^Yiq-j^Rqi), 



(25) 
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2s = ^(s l l t "Lb j )^2e q {q :j -f^Rq l ), 
i 

Qio 



16tt 2 

16tt 2 

e 2 
16tt 2 

„2 



m,bS l (j tlv RbiF^, 
{s^Lb){h»l), 



16n- 



■(s^Lb)(h^ 5 l). 



(26) 



To match the operator set in (16) onto these operators, we use the equations of motion to reduce all remaining 

LR and °Ib. 



two-quark operators to the gluon and photon magnetic moment operators 0\ R and 0\ R which are redefined as Q&g 



and Q 77 in (26). 





a) (b) (c) 

FIG. 3. Charged boson mediated penguin diagram contributing to b — > S7. 

Firstly, we present the values of the Wilson coefficients at raw scale calculated in NDR and HV schemes as follows 
with small V u bV* s omitted: 

Ci(m w ) = —(1 + 20 — , 



C2{m w ) = 1 



~0b 



K tb Kl 6 



11 (l+2() a ^ mff ' 35 a 



4tt 



18 4tt 



C 3 (mw) 



1 z. 



1 a s (m w ) 

QK tb K* s 6 4tt 

C^m^) = -3(75(771^) = C 6 (tovv) = 

2 . , 



E{x t ) + -£ 



q em 2£(x t ) +C(a; t ) 
67T sin 2 1 



E(x t ) + -t 



a s (m w ) 
4^ : 



C 7 (m w ) 



■ sin »n' 



Cs^w) = Cio(m w -) = 0, 



C 9 (m w ) = -^ cos2 #w 



C 7l {m w ) = --A(x t ) + 



Zsb 






67T 


= o, 




Zsb 




K tb K* ts 


6tt 


Zsb 


'23 


K tb K ts * 


36 ~ 



iC(x,) + D(x t ) --(! + () 



sin 6V 

4) 



9 



-f,,i , !///,, i +- C 3 j? (to+ ) - ic p i.2(m^ ) + ^C f m 



1 1 • 2„ 

-< ,, [ - + - sin W 



(27) 
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C 8G {rn w ) = - 1 -F{x t ) + -^— 
2 K tb K ts * 



C 9 (m w ) 
C w {m w ) 



-(-1 + 4sin 2 



OLera KtbK? a 

TT z sh B(x t ) - C(x t ) 



G{ Xt )-B{x t ) _ AC{xt) _ D{xt) + l {1 + + O 



Otem K t bKt s 



O 



sm &w 



sin tiw 

K tb K* tt 



(28) 



where £ = in NDR scheme and £ = — 1 in HV scheme. The loop functions can be found in Appendix B. 




FIG. 4. Tree level Feynman diagram contributing to b — > S7 and b — > si + Z . 

The terms proportional to K **x* in £^7,9(771^/) and the first term of Cg,io(mw) come from the tree-level diagram 
as displayed in Fig. 4, the term proportional to k*^k* m ^(toiv) comes from tree diagram due to the non-unitarity 
of CKM matrix in VQM. In expressions of Ct^sg^w); the first constant terms proportional to K **x* come from the 
charged current one-loop diagram Fig. 3 due to the non-unitarity of CKM matrix in VQM, other terms proportional 
to XtbK* come from the neutral current one-loop diagrams shown in Fig. 5. As in expressions Cg i io(mn/), the second 

terms proportional to denoted as O ( k **k* j come from the charged and neutral current one-loop diagrams. 

In deriving above equation, we have used the unitarity relation 

ZibZ* As = -z sb (29) 

which is a direct result of Eq. (4). 



S S 




FIG. 5. Neutral boson mediated penguin diagram contributing to b — *• sj. 
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At this moment, we would like to point out that the contributions from the Z-penguin charged current one-loop 
diagrams in the VQM to Cg^o(mw), have divergent terms due to the non-unitarity of CKM matrix. Although these 
divergences can be removed by renormalizing the tree level FCNC which exists in VQM Lagrangian, they are scheme 
dependent. However, these scheme dependences are subleading compared to the first terms of Cg t io(m\v), which we 
will neglect them in our calculation. Under this approximation they are not relevant for the inclusive dilcptonic 
decays. 



B. Scenario B: 8 D = Z ^^ H ~ 1 



In this scenario, the top and D quark, W and Z bosons can be integrated out together. The corresponding initial 
values of Wilson coefficients CV 7 ,8G( m w) are changed to [10] 



1 



C 7l (m w ) = --A(x t ) + 

2 KtbJ<ts 



1 



Zsb 



C 8G (mw) = —zF(x t ) + 

2 K tbJ< ts 



23 
36 
1 
3 



+ (fE(VD) + f§ (wd) + fUvo)) + e d ( i + i sin 2 0„ 



- 3 (fE(y D ) + f§ (wd) + fUvo)) ~ 3e d ( - + - sin 2 9 W 



(30) 



where yry = m 2 D /m 2 z , wr> = m 2 D /m 2 H . Other Wilson coefficients are the same as Scenario A we discussed. The 
functions stands for the contribution from boson x mediated penguin one-loop diagram with quark y in loops and 
are presented in Appendix B. 
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FIG. 6. The down-type vector quark contributions to Wilson coefficients of Qj 1 and Qsg at mw scale in unit of K 
(m,H — m t ). The solid and dashed lines stand for the results in Scenario A with or without consideration of QCD running 
effect from mc to mw, respectively. The corresponding values in Scenario B are denoted by dot-dashed lines. 



In Fig. 6 we plot the down-type vector quark contribution to the Wilson coefficients CV 7 and Csg as a function mo 
at mw scale, and demonstrate the QCD running effects from mp down to mw in Scenario A. 

As a consistency check, we can see that if the QCD correction is neglected by setting £ = 1 and only leading 
order terms of 8d are kept, for on-shell quarks and photon, our results C 7l (mw) and Cg,G( m w) m Scenario A would 
produce the results exactly in Scenario B. 

At the end of this subsection, we make some emphasis on the results obtained by using different approaches as 
follows: 

The effect of running from mo down to weak scale on the contributions from neutral current mediated by D quark 
in Scenario A is large, as shown in Fig. 6. However, at mw scale, since dominant new contribution comes from the 







charged current diagrams, the total results of CV 7i 8g(w*w) are changed slightly. This indicates that the dependence 
of the branching ratio of B — > X s j on mjj is quite weak in both Scenarios. Therefore, we can not extract the mass 
of down-type vector quark from this analysis. 



IV. CONSTRAINTS ON Z SB IN VQM FROM INCLUSIVE DECAYS B ->■ X s i + r 

A. Solutions of Wilson coefficients 



To obtain some predictions for inclusive B rare decays, wc need to determine the Wilson coefficients at rrib scale. 
Expanding anomalous dimension matrix and Wilson coefficient as 



7 = 7 



4ir 



4ir 



) 2 + 



c( M ) = c^V) + c (1 V)^^ + 

47T 



(31) 
(32) 



we solve the Wilson coefficients up to order a s (mw) according to the RGE (23). The relevant 10 x 10 one-loop 
anomalous dimension matrix among Qi_io, 7*-°^, and two-loop anomalous dimension matrices among Qi_io with 
277, 8G ano - ^9 needed in our calculations are collected in Appendix A. 

Using the 10 x 10 matrix in Appendix, and initial values presented in Section II, we obtain the solutions 



10 10 

i,i=l ^tb^ts , 



1-2 2 
Vi2 + -V u + - sin 2 6 w Vn + - sin 2 6 w V l9 
6 3 3 



(33) 



where r\ = a s (m\y)/a s (p), C^\l) = cf\mw)- The eigenvalues a; are obtained by diagonalizing the 10 x 10 



anomalous dimension matrix ^J-. 

2po 



v(0) 



2po 



At leading order, the solution of C 7l (rj) is given by 



C 7l (rj) = r]^C 77 (m w ) + ^ (rf$t - v^) C 8G (m w ) 

1GG ~ 777 V 



7G 7 



7GG - 777 

Pji 7G7 0Gi 



I OR ' U n, 7GG 



+VJ 1 



1 



2/3 2/3 2/3 a, - 



2/3 

277 

a/ 



2/3 



<i fc cf a). 



For Cg°' (77) , the corresponding solution is 



(34) 



4%) 



2 S 6 



10 "1 10 



where 



10 (0) 

i=l ^° 

10 (0) 

E'1,11 T>-1 



1 - 2 - 2 

V2 + 77 Vi 3 + - sin 2 w V l7 + - sin 2 w V l9 
3 3 



(35) 



(36) 



To obtain the scheme independent one-loop matrix element of O9 in VQM, calculation up to next-to-leading order 
(NLO) is necessary, as the case of SM. To obtain Cg(/x) , we frist present the NLO calculations for Ci(fi) as 
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C? ) (v) = V u 1 r 1 a 'V lk Gf\l) + V- 1 



v- — v- 1 

2/3o 



n ai — 1 



T %cf (i) 



a,- - a ; 



Po 



The solution of Cg(rj) in NDR scheme is 

Zsb 



10 ] 10 

+ 4sin2 ^ + ^ g ft( ^ - 1} J + M^y » - 1} 



K tb K* ts 

+ rC( E ^B( E 0_ 4 
sin 0w y 
+ E [W" - !) H + (Si + ffl£(x t )) (r/ a ' +1 - 1)] . 



(37) 



(38) 



The numerical results for parameters ai,pi,pi, si, qi arc follows: 



ai 


= (- 


-0.8994, -0.5217, -0.4230, 0.1457, 0.2609, 0.4087, 0.1304, 


-1.4034) 


Pi 


= (0.1648, 0.2424, 0.1384, -0.0073, -0.3941, 0.0433, 0, 0), 




Pi 


= (- 


-0.0248, 1.0004, 0.0853, 0.0207, 1.6602, -0.0690, 0.5223, - 


1.6974), 


si 


= (- 


-0.3554, -0.3579, -0.3617, 0.0072, -0.2009, 0.0490, 0, 0), 




m 


= (- 


-0.2701, 0, 0.0918, , 0.0059, 0, 0.0318, 0, 0), 




r N DR 


= (- 


-0.0292, -0.1960, 0.1328, -0.1858, 0.8997, -0.2011, 0, 0). 





(39) 

One can check that in case of z sb = 0, the results are the same as those in SM [26]. 

Note four-quark operators Qi_io can contribute to one-loop matrix clement of Og. Defining the effective coefficient 

Cc\ clS 



10 



c e 9 ff < si+r\o 9 \b >= c 9 ( v ) < si+r\o 9 \b > + E C 1%) < si + r\Qi\b >, 



i=l 



we obtain Cg ' {rj) in VQM as 



sieff _ Zsb 

9 ~ K tb K* ts 



— (-1 + 4sin 2 6 W ) + — Y j p l (v a ' +1 - 1) 
a em a(m w ) 



a s (m w ) 

+ 1) r? DR + (si + qiE(x t )) ( V a ' +1 l)] 



sin Uw 9 



Cf + 3Ci 0) ] + ; hci 0) + <7<°> + 3Cf + c<°> 



where s is dilepton mass squared normalized by m^, and 



-,c(0) 



3c[ 0) + a 



(0) 



2 ~<~ JL/ 3 "+" W "+" T ' 

4Cf> +4C| 0) +3Cf +C ( (0) 1 



V(0) , ,=V(0) 



T H0) = _ : 



The function /i(z, s) which includes light quark-antiquark pair contributions has an expression as 

,/ x 8, m b 8, 8 4 2,„ ,i 

M z, a ) = --ln---lnz + ---s--(2 + s)|l-^ 



In 



1-k-I 
2 arctan 



i7r for x = 4^- < 1, 
== for a; = 4— > 1. 

a; — 1 '' s 



(40) 



(41) 



(42) 



(43) 
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Now we consider the long distance contributions to Ct/^ . Apart from J/ty family, resonance 4> should be added to 
long distance part of Cg due to tree level Zsb coupling. Similar to previous analysis [27], the non-perturbative 
contributions can be parameterized as 



Y Res (s) = -J 



j,s(0). 



yeff 



T[</> -> £+fr}m^ 

C e 9 ff + Y Res ( S ). 

=(0), 



i=l 



r[^(is) -> £+r>j(!f) 



Here k = 2.3 is a phenomenological factor, |Tg luJ K| can be fixed from the J/5 , ,5'' data. Tg^ and Tg {u> correspond 

to the coefficients of h(^§,s) and h(^,s) in expression of Cg^ respectively, Tg^ = T b (°). Since experimental 

measurement for B — > X s <^> is not available so far, we set the phenomenological factor for <f) to be unity in our 
numerical calculation. 



(44) 



=(o) 



B. Constraints on z sb from inclusive decays B — ► -X" s ^ + i? 



With all Wilson coefficients at b quark mass scale ready, the invariant dilepton distribution for B — > X s £+^ can 
be expressed in terms of the effective Wilson coefficients defined above as 



H^^f'^ (1 - sfR (45) 

V 47T / 487T 3 



where 



Ro = 4 (l + \C e 7 f f \ 2 + (1 + 2a) (|C 9 e// | 2 + |do| 2 ) + 12Re(C% f C* ff *), (46) 

where is determined by Eq. (34) using the vectors /3 7 and /?g calculated in HV scheme. 

Now we constrain the parameter z sb from B — > A" s Z + /~. The branching ratio of _B — > depends on the 

tree level FCNC coupling z s & and \K tb K ts \. \K t t,K ts \ can be written in terms of \K cb K cs \ and z s b- In fact, from the 
quadrangle in bs sector, Fig. 1, one can derive the relation, 

\K tb K* ts \ = \K cb K*J + \K ub K* us \ cos^3 + \z s b\ costf, 

We use Eq. (47) to determine l-fT^-K^I in Eq. (45). Then, from the branching ratios of B — > one can derive 

the constraints on \z a b\ and 9. In the relation Eq. (47), for \K cb K cs \, we have used the experimental values shown 
in Table I. \K ub K us \ in Eq. (47) was neglected in the numerical calculation. Integrating Eq.(45), we exclude the 
resonances J/^, V E ,/ contributions by using the same cuts as experiments [8] 0.20 < |mj+j- — mj/^\ < 0.35, 0.15 < 
|m ; + ; - — mxs,i\ < 0.30. The corresponding 2a experimental bounds on the size of z sb and phase 9 from experimental 
measurement for B — > X s l + l~ arc displayed in Fig. 7. From this figure, we obtain 

M < 1-40 x 10~ 3 (95% C.L.). (48) 

One can see clearly that the absolute value of the allowed phase 9 becomes smaller as the size of z sb increases. 



TABLE I. Input parameters are used in calculation. All masses in unit of GeV. 



a s (m z ) 


0.119 




1/133 


m z 


91.19 


A 2 


0.12 GeV 2 


mw 


80.41 


Br ex (B-> X c ev) 


0.104 


m p b oU 


4.8 


\K US \ 


0.2196 ±0.0026 


mf le 


173.8 ± 5 


\K CS \ 


0.97 ±0.11 


m c /m b 


0.29 ±0.2 


\K ub \ 


(3.6 ±0.7) x 10" 3 


m s /m b 
sin 2 9w 


0.02 


\K*\ 


(4.12 ±0.2) x 10~ 2 
1.166 x 10~ 5 GeV- 2 


0.231 


G F 
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4 



1a bound 
2o bound 



2 



CD 



-1 



■2 



■3 




-4 



0.02 0.04 



0.06 



0.08 



0.1 



0.12 0.14 
x 10 "' 



FIG. 7. The (|z s i,|, (?) contour in Scenario A constrained by Br ex (B — > X s l + l ). The dashed, dotted lines correspond to 
la and 2a experimental bounds of B — » X s l + l~ in (1), respectively. The solid line denotes the experimental central value 
6.2 x 10" 6 . The region between dot lines is allowed at la level. 



V. SOME PREDICTIONS FOR B -> X s £ + r AND B X s j 

In this section, subject to the constraints on z s \, from B — » X s / + Z _ , we will present predictions for the invariant 
dilcpton mass distribution, FB asymmetry, the zero point of FB asymmetry of B — > X s l + l~ , and the branching ratio 
of-B^X^. 




7.5 2 2.5 3 

m„(GeV) 

FIG. 8. The differential rate of B — * X s l + l~ normalized by Br(B — > X c lv) as a function of dilepton invariant mass. 
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Firstly, we consider the invariant dilepton mass distribution of B — » X s l + l~ . According to the study in subsection 
IV A, it is expected to be sensitive to the parameter z s b- This feature is shown in Fig. 8. 

In addition to the differential and total branching ratio, the forward-backward (FB) asymmetry can provide crucial 
information on new physics. The normalized FB asymmetry distribution is defined as 



A 



I d6 Lease ' ^njcose) _ 3 

13 — j^ ggeggf) — - ~R~ Re[{sa 



eff 
9 



dsd cos 8 

Thus, the zero of the FB asymmetry sq in VQM is determined by equation 



Re 



(s C, 



eff 



c 



10 



= 0. 



2CV 7 )C* 



(49) 



(50) 




-0.5 



0.1 



0.2 



0.3 



OA 



0.5 

s 



0.6 



0.7 



0.8 



0.9 



FIG. 9. The FB asymmetry of B -> X s l + r distribution. 



It is very interesting to analyze how the zero of the FB asymmetry (so) is modified in VQM. Unlike to the case of 
SM where Cio is real, the coefficient Ciq is complex generally in VQM. Furthermore, the contributions to Cg^io from 
tree level FCNC diagram 



K tb Ki 



» \Ci ree (m w )\ 



Zsb 



(l-4sin 2 6l 



(51) 



indicate that Cio can have large imaginary part. Therefore, sq in VQM will have large deviation from that in SM. 

As an illustration, we plot the FB asymmetry as a function of s in Fig. 9 corresponding to different size of z s b- To 
show how the zero point sq sensitive to both the size of z s b and phase 9, we display the dependence of zero point of 
FB asymmetry on z s b in Fig. 10. Figs. 9 and 10 indicate that, subject to the experimental measurement for branching 
ratio of B — > X s l + l~, FB asymmetry distribution and the zero point of FB asymmetry are very sensitive to the 
parameters z s b and phase 9, especially in the region 0.6 x 10~ 3 < \z s b\ < 1.2 x 10~ 3 . 
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0.02 



0.04 



0.06 0.08 



0.1 



0.12 
x 10 



FIG. 10. The (\z s b\,so) contour in VQM subject to the la bounds of Br ex (B — > X s l + l ). For specified \z 8 b\, the phase 8 
effect on so is also shown. 



Now we study the radiative decay B — >■ X s ^. The Branching ratio of B — ■> X s j can be expressed as 



Br(B - X s7 ) = Br^B -> X c ez> e )i?, uarfc ((5 r ' 



1-6 1- 



A 2 



(52) 



where z = ml/ml, g(z) — 1 — 8z + 8z 3 — z A — 12z 2 lnz is the phase space factor. 6 max — 0.99 [5]. The A 2 - 
dependent term comes from the nonperturbative corrections to the semileptonic and radiative B meson decay rates, 
and A2 = {m 2 B , — m 2 3 )/4 ~ 0.12 GeV 2 . The ratio R qua rk is defined as 



T[b -» x sl ] E -<>^- 5 ^/ 2 6 



K, 



rb 



F(\D\ 2 +A), 



T[b —> X c eu e ] ng{z) 

where function F is related to the next-to-leading QCD corrections to the semileptonic decay 

n -1 



F= 1- 



8 a s (m b ) 



1 - 



2 a s (m b ) h(z) 



7T 



(53) 



(54) 



The |D| 2 term in (53) stands for the contribution of b — ► S7 while A term which is cutoff (5 dependent, includes the 
virtual and Bremsstrahlung correction to b — > -X" s 7 [5]. They can be written in terms of Wilson coefficients as 



r< +7V In — 



(55) 



i<j 



where 7^ = /3 7 ,; for j = 1 — 6, 777'' = 7 77 and 7g7 ; = 7G7- The explicit expressions of 7^, and fij(5) can be found 
Ref. [5]. 



*,(o) 



;,(0) _ 
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3.8 



1 a experment bounds 



3.6 



3.4 



3.2 



QQ 



,«.:.sJ.=.:.S<..- 



2.8 



2.6 



10 



12 



14 



Br(B^xrr) (ia 6 ) 



FIG. 11. The correlation of the branching ratios of B — * X s ^ and B — * X s l l~ predicted in VQM. The solid, 
dashed, dot and dot-dashed lines in the figure stand for different \z 8 b\- Each \z s b\ corresponds to one of the five specified 
3 7r). The experimental la bounds denoted as rectangular box are taken from (1) and (56). 



Fig. 11 shows correlation of the branching ratio of B — > X s ^ with that of B — > X s l + l~ in VQM. It is obvious from 
the figure that Br(B — > X s 7) is not so sensitive to the phase, which is not the case for Br(B — > Within the 

experimental bounds of £? — > A S Z + Z~ , corresponding branching ratio of £? — > A s 7 predicted in VQM is consistent 
with the experiment measurement [29] 



Br ex (B -> X sl ) = (3.34 ± 0.38) x 10" 



(56) 



VI. B —> Kit DECAYS 



Having established constraints on b — > s FCNC from b — ► sll and b — ► 57, we evaluate _B — > A"7r decay rates. The 
large electroweak penguin contribution to b — > s transition has been suggested in the present data of B factories. 
Because in VQM, the tree level FCNC may generate the electroweak penguin contribution, the purpose of this 
section is that with FCNC constrained by rare decays b — > sl + l~, we study whether the large electroweak penguin 
contribution in B — > Kir can be explained or not. The correlation between di-leptonic decay and non-leptonic decays 
is characteristic feature of the present model. 

We start with the following effective Hamiltonian: 

Heft = ^ j {CI Qi + C 2 * Ca) + tin {CI Q\ + C* 2 Q% ) - 6 £ (C* &) j • (57) 

The operators Qi is the charge conjugation of Qi defined in (26), = K is K* bl and 

6? = (fiiT^Xu^Lsi), e^^Lui)^,,^). (58) 

Here we must keep the £ u term because there are tree level contribution to B — » A"7r decays. As seen from the 
effective Hamiltonian, tree-level Z FCNC generates the electroweak penguin contribution which is not suppressed by 
a em . In terms of isospin amplitudes, the electroweak penguin operators contribute to both AI = and AI = 1 
components. Because for AI — amplitude, the large strong penguin contribution is expected, new physics may 
dominantly contribute to AI — 1 amplitudes. 

Below, we estimate the B — ► Kit amplitudes within the factorization approximation. The four quark operator with 
the flavor structure b — » sss and b — ► sec do not contribute to the processes. We obtain 
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A(B+ -+ K\+) = iG F {& (04 - + (2a 6 - a,)R K ) M B ^ K 

+ [C«fl2 - (,t (04 + 010 + 2(a 6 + a 8 )i? B )] M K *' B } , 

A(B+ - X+tt ) = 



di 3 . . 



V2M BK ' n 



A(B° 



+ [a 2 £u - it (04 + 010 + 2(a 6 + a 8 )i? K )] M Bv ' K 
- [C«a 2 ~ 6 (04 + 010 + 2(a 6 + a 8 )R B )] M K ^ B ) , 
= iG F { [i u a 2 - i t (04 + aio + 2(a 6 + a 8 )i? K )] M Bv,K 
+it la i -^f + (2a 6 -a s )R B )M K ^ B }, 



(59) 

(60) 
(61) 



A(S° -» *V) = z^§ { 



di 3 . . 
Uy - &^(-«7 + ag) 



V2M 



+it (a i -^ + (2a 6 -a 8 )R K )M B ^ K 
-it (a 4 -^f + (2a 6 -a s )R B )M K ^ B } 



(62) 



where a 2t ^ = C* 2l _ Y + C* 2l /N, a 2t = C* 2l + C* 2t _jN. And 



R K = 



m 



K 



(m b - m u )(m s + m u ) ' 
M ab ' c = f c (M*-m 2 b )F^ b (m 2 c ). 

In derving the amplitudes, we have used 

<n°(p'))\(bu) v \B + (p) > = (p + pX 
<K+(p'))\(bs) v \B+(p) > = \(p + p% 
<K+(-p)ir°(p'))\(us) v \0> 
where q = p — p' , and 



R L 



(m b + m u ){m s - m u ) " 



(63) 



(P + J>% 



m% 




ml 









m% 




m\ 




q2 




m\ 




ml 



q t , 
■q„ 
q fl 



F B 



V) 



^q,Ft-"(q% 



q z 



jrl< .: 



(q 2 ) 



^q»F«^(q% 



(64) 



< 0\(su) A \K+( Pk ) >= if K Pk^, < 0\(bu) A \B(p B ) >= ifsPBu, 

< K + (p k )\(su) A \0 >= -if K Pkn, < 7r°(p 7r )|(uu - dd) A \0 >= -iV^f^p^. 

It is easy to see that the isopsin relation 

V2A(B+ -» K+tt°) + A(B+ -» K°tt+) = V2A(B° -» K\°) + A(B° -» K+tt'). 



(65) 



(66) 



Since there are many uncertainties in calculating B — ► .KV such as strong phases and formfactors, instead of the 
branching ratios of B — > i^7r, it is more reasonable to use their ratios to constrain the new physics effect. In this 
work, we use R c , R n and their difference to constrain Z FCNC. R c and R n arc defined as 



Rc 



2Br{B+ -> if+rr ) + 2Br(B~ -» if-jr ) 
Br(B+ -> K°tt+) + Br(B- -> K°Tr-) ' 
Br^ -» isT+TT-) + Br(B° -> 
2Br(B° ifOrr ) + 2Br(B° -> i?°7r°) ' 



(67) 
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In our numerical calculations, we take large N limit and use f„ = 132 GeV. Since q 2 = rn 2 K ,rn 2 . is rather close to 

the point q 2 — 0, we neglect the q 2 dependence in the formfactors, while F^~* K = 0.33, Fq^ 71 = 0.379/\/2- The 
SU(3)-breaking effects in the formfactors are also neglected. The contribution from the terms proportional to fg is 
also omitted because it works as a suppression factor fs/Mg [16]. 

Running quark masses appear in the matrix elements of (S — P)(S + P) penguin operators through the use of the 
motion equation. The running quark masses at p ~ m b = 4.8 GeV scale are given by 

m b (m b ) = 4.34GeV 7 m s (m b ) = 0.09GeV, m u (m b ) = m d (m b ) = 0.03GeF. (68) 

The NLO Wilson coefficients in NDR scheme at m b scale can be obtained from (33) and (37). 

As discussed in Section II, the phases 63 and are defined in (47) and are shown in Fig. 1. We can write #3 as, 

3 ~ -atan f^J , (69) 

where r\ and p can be detremined from the measuerments of CP violation and mixings of K and Bd systems. By 
assuming that the effect of tree level FCNC on K and Bd is small, we can use the standard model values for ij and 
p. The allowed region of #3 is: 

100° < 3 < 140°. (70) 

If the tree level FCNC contribute to K and Bd system, the allowed region of 03 can be significantly changed from 
Eq.(70). Therefore, in our numerical calculation, we relax this condition and scan the region for 3 from 0° to 360°. 
Now we can see the physics quantities are deterimined by three inputs: the size of Z FCNC coupling \z sb \ 1 the phases 
of 03 and 0. We plot the i? n and R c correlation in Fig. 12 with three specified values of 03. For given z sbl the 
allowed region changes with 03. For \z sb \ < 1.5 x 10 -3 and 03 < 45°, no region is allowed by R c and R n experimental 
measurements at 2a level. When 03 = 100° ~ 120°, as shown in the figure, there are larger allowed regions. 

To see this more clearly, in Fig. 13 we show the dependence of R c — R n on the phase 03 with typical values 
\z sb \ = (1.0, 1.5) x 10~ 3 and = (— 7r, — ^, 0, §,7r). The experimental measurements [30] 

R c = 1.15 ±0.12, Rn = 0.78 ± 0.10, R c - R n = 0.37 ± 0.16 (71) 

at la and 2a level are also displayed. 
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FIG. 13. R c — Rn as a funition of the phase #3. 
From this figure, we note that R c — R n is sensitive to the phase #3 and 9. 

The constraints on the contour (|z s f,|,0) from R c — R n are shown in Fig. 14. One can see that, for 6*3 = 110°, la 
experimental bounds restrict \z s b\ and 9 as \z s b\ > 0.9 x 10 -3 , \9\ > 1. Corresponding to 2a experimental bounds, 
smaller \z s b\ is allowed as \z s b\ > 0.2 x 10~ 3 . 



CD 



Q 3 =11(f 



R c -R n 1 a bound 
R c -R n 2 a bound 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

I* J 

FIG. 14. The (|z s ;,|,0) contour constrained by R c — R n at la and 2a level. 



Finally, we obtain the allowed regions from R c — R n and B — » X s l + l~ by varying 9^ as a free parameter. The 
bounds are shown in Fig. 15. At la level, there are no overlapping region allowd from B — > X s l + l~ and B — > Ktt. 
From Fig. 15, we obtain 
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0.2 x 10~ 3 < \z sb \ < 1.2 x 10" 3 , (95% C. L.) 



(72) 



where large phase \9\ is favored. Now we would like to point out that if #3 is in the region of [100°, 140°], only the 
region with negative phase 9 in Fig. 15 is allowed. 



B-R n 2 a bound 

c n 
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0.1 0.12 0.14 

x 10 



-2 



FIG. 15. The (|z s i,|,#) contour constrained by combination of B — > X s l + l~ and B — > Kir at la and 2a level. The solid and 
dotted lines correspond to la and 2a bounds from R c — R n , respectively. The dash-dotted line and dashed lines correspond 
to 2a and la bounds from B — > X s l + l~ . The allowed region with > corresponds to 180° < #3 < 360° whereas the allowed 
region with 9 < corresponds to 0° < #3 < 180°. 



VII. CONCLUSIONS 



In this work, we have studied the new physics effects on the B meson decays B —> X s ^/ 1 B — > X s £ + £~ and B — > Kir 
in the Vectorlike Quark Model. The QCD running effect from the mass of the down-type vector quark D to weak 
scale has been taken into account. We extracted rather stringent constraints on the size of the tree level Z FCNC 
coupling and its CP violating phase from B — > X s £ + £~ . Within the bounds, we investigated various observables such 
as FB asymmetry of b — > sl + l~, the decay rates of B — > X s ^ and B — > if7r. We found that the zero crossing point 
of FB asymmetry so of £? — > is very sensitive to both the size and phase of z s f, and can be useful in probing 

the new physics, and both experimental measurements for B — > X s l + l~ and B — > Kir decays can be explained within 
the framework of vector-like quark model. The upper bound on the size of Z FCNC comes from B — > XJ + l~ , 
while the lower bound, from B — > Kir decays. Considering that the B factories such as BaBar and Belle are running, 
measurements for inclusive and exclusive B decays with high precision are expected. Therefore, the VQM will be 
tested in near future. 
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APPENDIX A 



In this appendix, we present all anomalous dimension matrices needed in solving the Wilson coefficients. 
(1) For the mixings of operators Q^r, Rl H , fl^i we obtain 



We also obtain the same matrix for Q\ R , R 



and P}' A , P'l, we obtain the following result: 




as Eq.(Al). For the mixings among the operators 1 LR (i 



7 = 



8tt 2 
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which is in agreement with Ref. [22]. 

(2) The 10 x 10 one-loop anomalous dimension matrix among <2i_ 
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is given by [25]: 
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where the color number N = 3 is used. For b — > sl + l , f = u + d, u = 2, d = 3. 7 ( -°- ) is scheme independent. 
In NDR scheme the 10 x 10 two-loop anomalous dimension matrix among the Qi-io is given by 
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7 , 646(M-d/2) 











_ 22/ 
3 9 
■225 



106(M-d/2) 
243 



10 3 6/ 

-ll(M-d/2) 
3 

106(«-rf/2) 
-1354 81 (u-d/2) 
826(m-< 3 /2) 



-65 
9 

-1192 

243 
-2384 71/ 



502/ 

243 
61/ 

U576/ 
243 



243 
-130 
9 

II _ 

3 

-225 



81 
1 + 



83(n-rf/2) 
9 

704Q-rf/2) 

243, , , 
1192 _ 71(u-d/2) 

243 9 
65 V _ 502(«-d/2) 



71 



-99 



+ 4/ + 



22/ 

3 

68/ 
9 



81 












+ 



81 



L _ 2/ 
9 

+ 2/ 
+ 3 












1+2/ 
2 ~ 3 
-21 2/ I 
2 9 / 



243 
\ 
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(3) The mixing entries among Qr 7 .sG an d Qi-io are follows: 

32 32 28 

777 = y 7G 7 = — jj") 7GG = y- 

In HV scheme, the entries of 01^q Q can be extracted from 0^~£ [23] by substituting 

1 — > ^e d , u^^e„u, d^^e d d, 

and thus, they are given by 

/ 416 464 .416 232 „ 32 ,448 200 16 ,448 , 100 232 ,416 , 116 , 2 

(4) The mixing entries among Cg and <2i_io are follows [25] with same substitution in (A6): 

(o) / 16 -16 -16, d 1, -16, d „, -16, d, 
7 i;11 = — , — (u ^ - 3). — (« - 2 - 3 )' — - 2 } ' 

16, d, ,16 2 4 j2 . ,16 2 4, 2 , ,16 2 4 j2 8, ,16 2 4 j2 8,\ 

(u ), -( — u 2 + -d 2 ), -( — u 2 + -d 2 ), -( — u 2 + -d 2 + -), -( — u 2 + -d 2 + -) , 

9 V 2 7 ' v 3 3 ; ' v 9 9 ; ' v 3 3 9 h v 9 9 3 1 J ' 

ndr,(\)_( 64 1600 3712 64(u-d/2) 64 /2240 512 \ 64(u-d/2) 3520 3392 

7l ' n ~ V Y' 243 '~ 243 3 ' T ~ V^J" + 243 J' 3 ' 243 U ~ 243 ' 

64 2 16, 2 . .3520 2 1696 2 1856 64 2 16* 32 2240 2 256 2 \ 

- { y u + Y d)A ^s u + w d) '!43-- ( T u + y d) '-y- ( ^r u -243 d } J' 

7S11 = -2/9b, 7 1 ^f' (1) - -2A, 
where A) = H - 2//3,/?i = 102 - 38//3. 



APPENDIX B 



Some one-loop functions needed in our calculations are follows. 

(1) For the calculation of the initial values of Cis and Cg ; io at mw scale: 



A(x 
B(x 

C(x 

D(x 

E(x 

F(x 



8x 3 + 5x 2 -7x x 2 {2-3x) 
12(1 -x t f + 2(1 -a;) 4 

X X 

+ — tttt ma;, 



In a;, 



4(.t - 1) 4{x - l) 2 
~ 2 - 6x 3x 2 + 2x 



(,■-1)2 

3a; 4 



(x-1) 



+ 



hix, 



30a; 3 + 54a; 2 



19a; 3 - 25a; 2 
" 36(x- l) 3 18(a;-l) 4 
2, x(18- 11a; -a; 2 ) 

'3 ln " + 12(1 -a; 3 ) ' 



32a; + 8 



In a;, 



5x 2 + 2x 



3x 2 



In a;. 



4(1 -a;) 3 2(1 -a;) 4 
(2) For the calculation of one-loop diagrams with down-type vector quark D: 
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fE(x) = 



5x 2 



5x - 4 x(2x - 1) , 

+ T~7 TTjlnX, 



72(x-l) 3 12(x-l) 4 



fS ( x ) = - e dX 

f%(x) = e d x 



7x 2 - 29x + 16 'Sx - 2 



48(x- l) 3 
5a; 2 - 19a; + 20 
48(a;- l) 3 



+ 



8(x- l) 4 
x — 2 



In a; 



(*-!)■ 



In a; 



(B2) 
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